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Abstract 

We present a systematic expansion of all constraint equations in canoni- 
cal quantum gravity up to the order of the inverse Planck mass squared. 
It is demonstrated that this method generates the conventional Feynman 
diagrammatic technique involving graviton loops and vertices. It also re- 
veals explicitly the back reaction effects of quantized matter and graviton 
vacuum polarization. This provides an explicit correspondence between the 
frameworks of canonical and covariant quantum gravity in the semiclassical 
limit. 
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1. Introduction 



Attempts to understand the quantum nature of the gravitational field typically 
fall into two classes. One possibility is to first construct a well-defined mathematical 
framework for the full theory, from which eventually all quantum gravitational phe- 
nomena can be derived. Superstring theory and canonical quantum gravity in the 
Ashtekar formulation belong to this class. The other alternative seeks to apply vari- 
ous approximation schemes to a set of formal quantum gravity equations. One hopes 
that in this way some of the profound conceptual problems of quantum gravity and 
quantum cosmology - such as the problem of time and the related problem of Hilbert 
space - can be tackled. One also hopes to derive in this way quantum gravitational 
corrections to standard physics. Although thus in a sense less ambitious than the first 
class of alternatives, it is far more straightforward to get definite results which even 
may provide the first window towards genuine tests of quantum gravity. 

One such approximation scheme is to perform an expansion of the Wheeler-DeWitt 
equation of canonical quantum gravity with respect to the inverse Planck mass squared, 
trip 2 . This expansion was pioneered in Jl|, and later developed from various points 
of view. A most remarkable result is the recovery - at order m P - of the (functional) 
Schrodinger equation on a classical background spacetime as an approximate equation 
from quantum gravity || |, ||. At this order, however, two important aspects of the 
full theory are not yet incorporated: Back reaction effects of quantum matter on the 
dynamical gravitational background and the proper quantum effects of the gravitational 
field itself. Their description requires higher-order iterations of the Wheeler-DeWitt 
equations in inverse powers of the Planck mass. Such higher-order correction terms 
were derived in various forms in both the canonical theory 0, || and in the covariant 
framework of path integrals |10, 11], |12 |. 

In ||, definite correction terms to the Schrodinger equation were derived at order 
nip 2 from the Wheeler-DeWitt equation in a formal way (without addressing the issue 
of regularization) . In particular cases it was possible to evaluate these correction terms 



to find quantum gravitationally induced energy shifts |9|, |13j. It was also shown that 
these correction terms can in principle lead to observable effects in the spectrum of the 
cosmological background radiation JTjj] . 

Apart form the issue of regularization, the analysis in H was incomplete for two 
reasons: First, it applied the expansion scheme only to the quantum Hamiltonian 
constraint, but not to the full, complete, set of constraint equations. A complete 
picture has to implement all constraints which are interconnected by their commutator 



algebra [15]. Second, although a general expression for all correction terms to the 



Schrodinger equation was given in M , the final expression of these terms only involved 
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the so-called "longitudinal" derivatives - the derivatives with respect to configuration 
space coordinates along the classical trajectory. Only this part of the correction terms 
follows from the previous order equations alone, while the remaining terms encode the 
dependence on boundary conditions. The main point of the present paper is to present 
a concise method of giving explicit expressions for all correction terms (including the 
"transversal derivatives") for the full set of quantum constraints. It is shown that 
this leads to the conventional Feynman diagrammatic technique involving graviton 
propagator, vertices and loops and thus provides a concrete physical interpretation of 
all terms. 

One of the widespread approaches to canonical quantum gravity is characterized 
by the fact that it does not contain a notion of spacetime at the fundamental level 



(see e.g. [[L6| for a conceptual discussion). Spacetime emerges as an approximate 
notion on the level of the semiclassical expansion described above. But even within 
such a conceptual framework it would be desirable on this level to possess a spacetime 
covariant description of all correction terms. This would yield a bridge between the 
canonical formalism and the covariant effective action formalism [1C, 0. It would be of 
crucial importance for the correct covariant regularization of the ultraviolet divergences 
which inevitably arise in all loop orders of the l/m|>-expansion. This regularization 
should maintain the general covariance in the form not splitted by a (3+l)-foliation. 
The general covariance is not manifest due to the canonical origin of the constraint 
equations, but it gets restored in the proposed calculation technique, because the arising 
Feynman diagrams can be cast into spacetime covariant form. Another important 
advantage of this technique is that it clearly reveals the back reaction effects of quantum 
matter and graviton vacuum polarization. This back reaction has the form of special 
nonlocal gravitational potentials (similar to the retarded potentials in electrodynamics) 
contributing to the kinetic and potential terms in the effective Hamiltonian of the 
corrected Schrodinger equation. 

Our paper is organized as follows. Section 2 gives a brief introduction into the 



formalism of the Planck mass expansion. In Section 3 we review the work in |1| [19 



presenting a consistent factor ordering of the theory. We then give the calculation of 
the one-loop gravitational prefactor. In Section 4 we present a detailed calculation of all 
correction terms to the matter Hamiltonian up to order rrip 2 . As we shall describe, there 
will be three types of correction terms: a purely quantum gravitational term, a term 
corresponding to the contribution of quantum matter, and a cross term. Sections 5-7 
give a description of all the correction terms in the language of Feynman diagrams 
(some computational details are relegated to appendices). Section 8 gives an account 
of the back reaction of matter on the gravitational background as well as a comparison 
with the terms calculated in ||. Section 9 contains a summary and an outlook on 
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future work. 



2. The Wheeler-DeWitt and Schrodinger equations 

Consider the Lagrangian action of Einstein gravity theory coupled to matter, 
S**\gap(x)M*)] = m "p I d 4 x v ^ 4 i?(^(x)) + 5 mat [^(x),^(x)] , (2.1) 



where rnp = l/16nG, and the cosmological constant is set equal to zero. The canonical 
form of this action in terms of three-metric coefficients g ab , matter fields if and their 
conjugate momenta p ab ^p^ (a, b, ... = 1, 2, 3) has the form |2(| |2l 



5 tot [^ fc ,p afc ,^,^,iV ± ,iV a ] = J dtd 3 x (m 2 P p ab g ab +p^ - iV x T ± - N a T a ) , (2.2) 

where the lapse N 1 - = (— ^y 00 ) -1 / 2 anc i shift N a = g abll gob functions enter as Lagrange 
multipliers of the gravitational constraints 

T ± = m 2 P H ± (g ab , <p, p ab , p v ) + (g ab , <p, p ab , Pip ), (2.3) 
T a = m P H a (g ab ,v,p ab ,p v ) + Hr t (gab,V,p ab ,Pv)- (2-4) 

The gravitational parts of these constraints are given by the gravitational superhamil- 
tonian and supermomentum 

H ± (g ab , v,P ab ,P v ) = \G ab , cdP ab p cd - Jg% (2.5) 
H a (g ab ,^p ab , Pv> ) = -2g ab V c p bc } (2.6) 

where G abjCd = g~ l/2 {g a c9bd + 9ad9bc - 9ab9cd) is a local DeWitt supermetric and V c 
denotes the covariant spatial derivative. The matter parts of the constraints, H™ at and 
H™ at , depend on the concrete choice of matter action which we shall not specify here. 
Its form can be strongly constrained from general principles such as ultralocality |22 . 
The above constraints 

T ± « 0, T a w 0, (2.7) 

are in Dirac's terminology first class constraints having the following closed Poisson 



algebra flll g|, gg 

{T ± (x),T ± (x')} = g*(x)T h (x)d tt 6{x,*!) - (x <-> x'), 
{T ± (x),T a (x')} = -T ± (x'R5(x',x), 

{T a (x),T 6 (x')} = T b (x)9 a (5(x,x / )-(a,x^6,x / ) . (2 A 



3 



These constraints imply that the phase space variables are not dynamically inde- 
pendent - a property which in the Dirac quantization framework leads to quantum 
constraints on physically admissible quantum states. In this quantization scheme the 
c-number constraints Q2.3| ) - ( |2.4| ) become the operators (T±,T a ) which should satisfy 
the closed commutator algebra generalizing ( [2.8| ) to the quantum domain. They select 
the physical quantum states | by the equations |], 

TJ*> = 0, r a |*) = 0. (2.9) 

In what follows we shall denote the quantum states in the standard representation 
(Hilbert) space of matter fields (<^(x),p ¥ ,(x)) by ket-vectors, the matter field operators 
label by hats and choose the functional coordinate representation for metric variables 

$«6(x) = <7a6(x), P" fe (x)= , \ \ (2.10) 

rnij, dg ab {-K) 

Note that we have redefined the gravitational momenta by a factor m? P compared with 
the standard convention. 

The quantum states then become the state vector- valued functionals of three-metric 
coefficients ^[^(x)]), and the equations ( |2.9|) take the form of well known Wheeler- 
DeWitt equations 

\-^G ab , cd —- ml^R + HT 1 ) mg ab \) = 0, (2.11) 

( 2mf> dg ab bg cd J 

"{ 2 5 ~ 1" 

-TfaVe— + H™ at mg ab }) = 0, (2.12) 
{ i og hc J 

where the inverted commas indicate that these functional variational operators are a 
symbolic representation of some operator realization of the classical constraints ( |2.5| ) - 
J), implying both the operator ordering and quantum corrections proportional to h. 
The form of and ( gjj) shows that the size of the quantum gravitational 



effects is governed by the dimensional parameter m 2 P - the square of the Planck mass. 
The quantum effects of matter fields (their quantum commutators, couplings, etc.) in 
turn are determined by U chosen for simplicity to be one in our units. The justification 
of this choice is that in our paper we shall develop the semiclassical expansion in 
1/mp of quantum gravitational effects without ^-expansion of quantum matter. The 
lowest-order approximation of such an expansion for solutions of the Wheeler-DeWitt 
equations |2], |3|, || yields the following semiclassical form 

m 9 a b ])=e [m2pS[gab] mg ab ]), (2.13) 
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where S[g a b] is a purely gravitational Hamilton- Jacobi function. This is a solution of the 
vacuum Einstein-Hamilton- Jacobi equations - the gravitational constraint equations 
with the momenta replaced by their Hamilton- Jacobi values (gradients of S[g ab \): 

H ± (g ab , 5S/5g ab ) = ]-G abiCd ^-^- - jfR = 0, (2.14) 
2 og ab bg cd 

H a (g ab , 5S/5g ab ) = -2g ab V c — = . (2.15) 

og bc 



Substitution of ( |2.13| ) into the Wheeler-DeWitt equations leads to new equations 
for the state vector of matter fields |$[g afe ]) parametrically depending on the spatial 
metric 

I -G ab c d- — h H™ at (g ab ) 

{ i dg ab og cd 

1 " 5 2 S " 1 " 5 2 "1 

+^G abtCd —- ——Gat.cd , , |$M) = 0, (2.16) 



2i og ab bg cd 2mf, og ab 6g, 



cd 



| - jfaVcj^ + Hr\g ab )} Mg ab ]) = o. (2.17) 

The conventional derivation of the Schrodinger equation from the Wheeler-DeWitt 
equations then consists in the assumption of small back reaction of quantum matter 
on the metric background which at least heuristically allows one to discard the third 
and the fourth terms in ( |2.16| ). Then one considers |$[<7a&]) on the solution of classical 
vacuum Einstein equations g ab (x.,t) corresponding to the Hamilton- Jacobi function 
S[g a b), 

\<S>(t)) = mg ab (x,t)\). (2.18) 

After a certain choice of lapse and shift functions (N- 1 , N a ), this solution satisfies the 
canonical equations 

g ab = N^Ga^d— + 2V (a iV 6 ), (2.19) 

so that the quantum state ( |2.18| ) satisfies the evolutionary equation obtained by using 

- Mt)) = I d 3 xg ab (x) -—-Jcl^]) (2.20) 
dt J og ab (x) 

together with the truncated version of equations ( |2.16|) - ( 2.17 ). 

The result in the approximation of the above type is the Schrodinger equation of 
quantized matter fields in the external classical gravitational field, 

i^|$(t)) = tf mat |$(t)>, (2.21) 

^mat = J d 3 x |]V ± (x)^ at (x) + iV" (x)fl^ at (x) } . (2.22) 
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Here, H mat is a matter field Hamiltonian in the Schrodinger picture, parametrically 
depending on (generally nonstatic) metric coefficients of the curved spacetime back- 
ground. Examples can be found e.g. in ||, |13| . 

Such a derivation of quantum field theory from the Wheeler-DeWitt equations 
dates back, at the level of cosmological models, to the pioneering work of DeWitt 
[IJ. It was later performed by Lapchinsky and Rubakov for generic gravitational 
systems and discussed in various contexts in ||, [|, ]5|, || [7], |], ||. This method turns 
out to be the most commonly used approach and, in fact, the most effective way 
of interpreting the cosmological wavefunction in the semiclassical approximation. It 
establishes in particular the links between fundamental quantum cosmology and the 
physics of the early Universe. On the basis of the obtained Schrodinger equation this 
method obviously generates all quantum effects of matter fields in curved spacetime, 
but in this order of approximation it does not yet contain two important ingredients 
of the complete quantum scheme: i) purely gravitational quantum effects and ii) back 
reaction effects of quantum matter on the gravitational fields. Obviously, these effects 
are generated by the discarded third and fourth terms of ( j2.16| ). The purpose of this 
paper is to develop the systematic method of 1/mp-expansion perturbatively generated 
by these terms and demonstrate that this method gives rise to back reaction effects 
and can be described by a Feynman diagrammatic technique including graviton and 
graviton-matter loops. 

From a methodological point of view, we shall show how the treatment of the grav- 
itational preexponential factor presented in |TT, [T2|, [18|, |1^] improves the results of [[7] 



where the purely gravitational effects were not properly disentangled. We shall also 
give transparent expression for the "transversal derivatives" of the correction terms 
||. As we shall see, these terms generate the kinetic part of the nonlocal gravita- 
tional potentials describing the back reaction of quantum matter on the gravitational 
background. 

There is an obvious problem arising in the implementation of this program. The 
effects we are going to consider go beyond the tree-level approximation and, there- 
fore, require the precise knowledge of the operator realization of the Wheeler-DeWitt 
equations ( |2.11| ) - ( |2.12| ) and ( |2.16| ) - ( |2.17| ), that is replacing the inverted commas 
by actual operator ordering of coordinates and momenta. The main property that 
these operators should satisfy is a closure of the quantum commutator algebra of con- 
straints generalizing ( |2.8| ) to the quantum level as consistency conditions of fl2.9p . In 
the coordinate representation this realization is known for generic theories subject to 
first-class constraints in the one-loop approximation Fortunately, for theories with 
constraints of the gravitational type quadratic and linear in momenta (like ( |2.5|) - ( [2.6|) ) 
this realization is known in the exact theory, formally closing the commutator algebra 
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beyond the semiclassical expansion |T9]]E} In the next section we dwell on the construc- 
tion of this realization in terms of the condensed DeWitt notation which will also allow 
us to formulate the necessary gravitational prefactor before going over to higher orders 
of the 1/mp-expansion for the Wheeler-DeWitt equations. 

3. Operator realization of the Wheeler-DeWitt equa- 
tions and the one-loop gravitational prefactor 

The operator realization of the Wheeler-DeWitt equations in the full theory not 
restricted by any minisuperspace reductions uses a very complicated formalism. What 
makes this formalism manageable is the use of condensed DeWitt notations which 
formally represents the complicated field system in terms of a quantum-mechanical 
model with a finite-dimensional phase space having a finite-dimensional space of local 
gauge (diffeomorphism) transformations. This can be achieved by introducing the 
collective notation for gravitational phase space variables 

q l = g ab (x), ft=p a6 (x), (3.1) 

in which the condensed index i = (ab, x) includes both discrete tensor indices and three- 
dimensional spatial coordinates x. Similar notations for the full set of constraints and 
their gravitational and matter parts yields 

T M ( g ,p) = (T L (x),T a (x)) > T,(q,p)=m 2 P H, + H™\ (3.2) 

H,(q,p) = (# ± (x), J7 a (x)), H™\q,<p,p ¥ ) = (flr*(x),fl?*(x)) ■ (3-3) 

The index \x enumerates the superhamiltonian and supermomenta of the theory, as 
well as their spatial coordinates \i — > (\i, x). In these notations the functional depen- 
dence on phase space variables is represented in the form of functions of (q % ,Pi), and 
the contraction of condensed indices includes integration over x along with discrete 
summation. In condensed notations the gravitational part of the canonical action ( |2.2j ) 
has the simple form 

S[q,p,N] = Jdt (p t q l -N»H,(q,p)) , (3.4) 

with the superhamiltonian and supermomenta given by expressions which are quadratic 
and linear in the momenta, respectively, 

H± = \GfpiP k + V ± , H a = V l aPl . (3.5) 



1 The formal closure of the commutator algebras of |T^, |T^] holds up to possible quantum anomalies 
arising when regulating ill-defined products of operators taken at the same space point. The problem 
of these anomalies in field models in spacctimc dimensions higher than two is far from resolution and 
goes beyond the present paper. 
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Here the indices _L— > (J-,x) and a — > (a, x) are also condensed, G 1 ^ is the ultralocal 
three-point object containing the matrix of the DeWitt supermetric, V± denotes the 
potential term of the superhamiltonian (|2.5|) , and is the generator of the spatial 
diffeomorphism of q l (see below). The objects G 1 ^ and have the form of the following 
delta-function type kernels |TjJ 

Gf = G ahfid 5(xi,x fc )5(x_L,x fe ), i = (o6,Xi), k = (cd,x k ), _L=(_L,x_l), (3.6) 
V„ = -2g ( 6 V c )5(x a ,Xi), i = (be, Xi), a = (a,x a ). (3.7) 

Note that the object G l ± itself does not form the DeWitt supermetric because it 
contains two delta-functions. Only the functional contraction of Gj_ with the constant 
lapse function N x = 1 converts this object into the distinguished ultralocal metric on 
the functional space of three-metric coefficients, 

G ik = G * N ±^^ _ J d 3 x± Q ik = Gab cd 5 ( x . )Xfe ). (3. 8 ) 

The Poisson bracket algebra for the gravitational constraints in condensed notations 
can be written as 

{H„H u } = U« u H a (3.9) 
with structure functions U" = U^ u (q) that can be read off from ([2.8] ) . This algebra 



implies the gauge invariance of the action ([3.4]) under the transformations with local 
(arbitrary time and space dependent) parameters = ^(t). These transformations 
are canonical (ultralocal in time) for phase space variables, 

<V = {q\ H^f, 5 Pi = {pi, H^f, (3.10) 

and quasilocal (involving the time derivative of J 711 ) for Lagrange multipliers p4| , 

5N» U£ v N a F v . (3.11) 

Note that the transformations of phase space coordinates 

OH, 

dpi 

have as generators the vectors which are momentum independent for spacelike dif- 
feomorphisms fi = a (and, therefore, coincide with the coefficients of the momenta in 
the supermomentum constraints (13.70 ), but involve momenta for spacetime diffeomor- 
phisms normal to spatial slices, = G±Vk- 

With these condensed notations let us formulate the operator realization of grav- 
itational constraints H^(q,p) — > closing the commutator version of the Poisson 
bracket algebra ( |2.8| ) 

H„H U ] =iU^H x . (3.13) 



5q l = V^, V; = ^ (3.12) 
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As shown in ||19H , the fact that (|3.13| ) holds follows from the classical gravitational 



constraints (|3.5| ) by replacing the momenta pk with the functional covariant derivatives 
T>k/imp - covariant with respect to the Riemann connection based on the DeWitt 
supermetric ( |3.8| ) - and by adding a purely imaginary part (anti-Hermitian with respect 
to the L 2 inner product): the functional trace of structure functions, ihU^/2. With 
this definition of covariant derivatives it is understood that the superspace of three- 
metrics q is regarded as a functional differentiable manifold, and quantum states {^(q)^ 
are scalar densities of 1/2-weight. Thus the operator realization for the full constraints 
including the matter parts has the form |fL9|l p| 

f ± = -J^-GfViVk + m 2 P V ± + ~Ul v + HT\ (3.14) 
2mp 2 

f a = W a V l+ l -U: u + H^. (3.15) 

Imaginary parts of these operators are either formally divergent (beeing the coincidence 
limits of delta-function type kernels) or formally zero (as in ( |3.14|) because of vanishing 
structure functions components). We shall, however, not ascribe to them particular 
values characteristic of gravity theory and keep them of a general form, expecting that 
a rigorous operator regularization will exist that can consistently handle these infinities 



as well as corresponding quantum anomalies (see the discussion of this point in [18|, |19|). 

As shown in |12|, [19]], the DeWitt supermetric on superspace has as functional Killing 
vectors the generators of spatial diffeomorphisms V„, T> l V k + V k V l a — 0, T> % — G lh V k , 
and in view of its ultralocality, G lk ~ 5(xj,Xfc), the covariant derivative conserves not 
only the metric itself, but also the three-point object G± , T> m G %k = 0. Therefore, 
the kinetic terms of operators ( [3.14|) - (|3.15|) do not require additional prescriptions 



of operator ordering and, moreover, turn out to be Hermitian with respect to the 
(unphysical) L 2 inner product. (This does of course not hold for the antihermitean 
contributions to the potential term.) 

An important result of jIT[ [12], |18| , which we shall need below, is a closed construc- 
tion of the one-loop preexponential factor for a special two-point solution K (q, q') of 
the Wheeler-DeWitt equations represented in the form of the semiclassical ansatz 

,, im 2 P S(q,q') 

K(q,q') = P(q,q')e P W ' . (3.16) 



2 We assume that the operators of matter parts of constraints i?™ at independently satisfy the 
commutator algebra with the same operator structure functions as in the gravitational sector ( |3.9| ), 
while the commutator between and 7?™ at is identically zero |2^| . With respect to possible quantum 
anomalies this assumption is rather nontrivial, but this goes, as it has already been mentioned above, 
beyond the scope of this paper. 
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Here, S(q, q') is a particular solution of the Einstein-Hamilton- Jacobi (EH J) equations 
with respect to both arguments - the classical action calculated at the extremal of 
equations of motion, joining the points q and q' in superspace, 



H^q,dS/dq) = 0. 



(3.17) 



The one- loop (0(m° P ) part of the 1/mp-expansion) order of the preexponential 
factor P(q, q') here satisfies a set of quasi-continuity equations which follow from the 
Wheeler-DeWitt equations at one loop [[II], [18], 



^(v;p 2 ) 



p = dS I dq 



(3.18) 
(3.19) 



with the generators here evaluated at the Hamilton- Jacobi values of the canonical 
momenta. The solution of this equation found in [II], [13], H5| turns out to be a particular 



generalization of the Pauli-Van Vleck-Morette formula p5| - the determinant calculated 
on the subspace of nondegeneracy for the matrix 

d 2 S(q,q') 



'ik' 



dq 1 dq k ' 



(3.20) 



This matrix has the generators as zero-eigenvalue eigenvectors [II], [I2|, [OH- An 
invariant algorithm of calculating this determinant is equivalent to the Faddeev-Popov 
gauge- fixing procedure leading to the solution of ( |3.18| ). It consists in introducing a 
"gauge-breaking" term to the matrix (|3.20| ) . 



ik' 



Sik' + 



i C fJ.U<Pk' > 



(3.21) 



formed with the aid of gauge-fixing matrix c^ u and two sets of arbitrary covectors (of 
"gauge conditions" ) 0f and at the points q and q' respectively, satisfying invertibility 
conditions for "Faddeev-Popov operators" at these two points 



J» = tfVl J = detJ/f ^ 0, J£ = 0^Vf, f = detJ'^ ^ 0. 



(3.22) 



In terms of these objects the preexponential factor solving the continuity equations 
(|3.18|) is given by the following expression |ll], |12|, |18 | 



detF 



■ik' 



J J' det c 



1/2 



(3.23) 



which is independent of the introduced arbitrary elements of gauge-fixing procedure 
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4. Semiclassical ansatz for the two-point solution 



We shall now proceed to perform the semiclassical expansion for solutions to the 
Wheeler-DeWitt equations. Since we are interested in giving an interpretation in terms 
of Feynman diagrams, we shall not consider wave functionals as in |J , but two-point 
solutions ("propagators"). Due to the absence of an external time parameter in the 
full theory, such two-point functions play more the role of energy Green functions 
than ordinary propagators |27|. However, in the semiclassical limit, a background time 
parameter is easily available, with respect to which Feynman "propagators" can be 
formulated. 

Let us therefore look for a two-point solution of the Wheeler-DeWitt equations in 
the form of the ansatz 

K(q + ,q-) = P(q+,q-)e U(q+,q-), (4.1) 

where we denote (as agreed above) by a hat the operators acting in the Hilbert space 
of matter fieldsQ. Here, S(q+,q-) is the principal Hamilton function satisfying the 
gravitational EHJ ( |3.17 ) equations, and P{q +) qJ) is the preexponential factor ( |3.23| ). 
Substituting this ansatz into the system of the Wheeler-DeWitt equations and taking 
into account the EHJ equations and the continuity equations for P(q+, <?-), we get for 
the "evolution" operator U(q + ,q_) the equations 

iV k ± V k U = HT'U - -^p-^GTV m V n {P U), (4.2) 

Zrrip 

iV k a V k U = H^U, (4.3) 

where all the derivatives are understood as acting on the argument q + . 
Evaluating this operator at the classical extremal q + — > q(t), 

U(t) = U(q(t),q_), (4.4) 
where q(t) satisfies the canonical equations of motion corresponding to S(q+,q~), 

q l = A^V;, (4.5) 
one easily obtains the quasi-evolutionary equation 

i^U(t) = H cS U(t) (4.6) 



with the effective matter Hamiltonian 

1 

i P 



H cs = ^mat -N ± G^ n 'D m T> n [PU]p- 1 U- 1 . (4.7) 

2?TZ i 



3 We note that in M 1/D corresponds to P, while % corresponds to U . 
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(Recall that this is an integrated equation, cf. (2.20).) The first term on the right-hand 
side is the Hamiltonian of matter fields at the gravitational background of (q, N)- 
variables if mat ; 

The second term involves the operator U itself in a nonlinear way and contributes 
only at order nip 2 of the expansion. Thus, ( [4.6[ ) is not a true linear Schrodinger 
equation, but semiclassically it can be solved by iterations starting from the lowest 
order approximation 



Uo = T exp 



-i f t+ dtH mat 



(4.9) 



H cti = ^mat_ ^ 

Here, T denotes the Dyson chronological ordering of the usual unitary evolution op- 
erator acting in the Hilbert space of matter fields (yj,^). The Hamiltonian if mat = 
H mat ((p,p v ,q(t),N(t)) is an operator in the Schrodinger picture of these fields (tp, p^) 
parametrically depending on gravitational background variables (q(t), N(t)), i.e. eval- 
uated along a particular trajectory ( "spacetime" ) in configuration space. 

The Dyson T-exponent obviously explains the origin of the standard Feynman dia- 
grammatic technique in the matter field sector of the theory which arises in the course 



of the semiclassical expansion of ( |4.9| ). We shall show now that the gravitational part 
of this diagrammatic technique involving graviton loops naturally arises as a result of 
iterational solution of ( |4.6| ) - ( [4.7| ) in powers of 1/mp. 

The effective Hamiltonian in the first order approximation of such an iterational 
technique can be obtained by substituting ( |4.9|) into (|4.7| ) to yield 



Hf(t+) = # mat - -^g™(v m v n u )u - 1 

2mp 

1 , 1 



2 g mn (v m v n p)p- 1 -—g mn (v m p)p~ 1 (v n u )u 1 . (4.n) 

2m P m, 



op flop 

Here we have used the new notation 



gmn = jylgmn 

for another metric on the configuration space (compare with ( |3.8|) ) which uses the ac- 
tual value of the lapse function corresponding to the classical extremal (f4.5|). This lapse 
function generally differs from unity. We also decomposed the first-order corrections in 
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the effective Hamiltonian into three terms corresponding to the contribution of quan- 
tum matter (generated by Uo), purely quantum gravitational contribution generated 
by P, and their cross termQ. 

Further evaluation of these terms demands the knowledge of derivatives acting on 
the configuration space argument q + of TJq and P. Obtaining these derivatives leads to 
the necessity of considering the special boundary value problem for classical equations 
of motion, the graviton propagator and vertices - elements of the gravitational Feynman 
diagrammatic technique. The rest of the paper will be mainly devoted to present the 
corresponding derivations. 

We begin by introducing the collective notation for the full set of Lagrangian grav- 
itational variables, which includes both the spatial metric as well as lapse and shift 
functions, 

0° = ( 9 *(t),W(t)) . (4.13) 

Actually this collection of ten fields g a ~ g a p{^, t) comprises the whole set of spacetime 
metric coefficients taken in a special parametrization adjusted to the (3+l)-splitting. 
In what follows we shall need also spacetime condensed DeWitt notations in which 
the index a includes not only the spatial coordinates x but also the time variable t, 
and contraction of these indices will imply also the integration over t. In spacetime 
condensed notations the gravitational action has the form 

S[g] = dtL(q(t),q(t),N"{t)) ■ ( 4 - 14 ) 

The Lagrangian does not involve the time derivatives of lapse and shift functions iV^t) 
which play the role of Lagrange multipliers for the first-class constraints. This La- 
grangian is related to the integrand of the canonical action 



L(q,q t N») = (ptf - N»H„(q,p))\ p=pW) (4.15) 

by the substitution of the expression for the canonical momentum p®(q, q, N) in terms 
of the velocities q - the solution of the canonical equation of motion 

■i = N n^±H . ( 4 . 16 ) 

OPi 

In the notations of the above type the evolution operator Q4.9p as a function of q + 
can be regarded as a functional of the classical extremal g a = {q l {t), N^it)) joining the 
points q + and g_ in superspace at the respective moments of time t + and 

U G {q + ,qJ) = U Q [g(t\q + ,qJ)\. (4.17) 



4r 



The second correction term in ( 4.11 ) was "absorbed" in H by the introduction of some wave 



functional a obeying a second-order WKB equation. 
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Here, g(t\q +) qJ) is a classical extremal of vacuum Einstein equations parametrically 
depending on end points q± at t±. In terms of the action ( [4.14 ) the boundary value 
problem for this extremal can be written as 

|^ = 0, (4.18) 
6g a (t) v ' 

q{t±) = q± • (4.19) 
Thus, the derivative of Uq with respect to q + can be written as 
dU ft* dg«(fy 5U [g] ^dg^5Uo 



dqf Jt- dqf 5g a (t) dq% 5g a 

Note that in the last equality we used the contraction rule for spacetime condensed 
indices, implying the time integration together with the spatial integration. In what 
follows we shall not introduce special labels to distinguish between condensed spatial 
and spacetime notations. As a rule, when the time argument is explicitly written, we 
shall imply that the corresponding condensed index or indices involve only spatial coor- 
dinates and their contraction does not involve implicit time integrals. Another distinc- 
tion between these two types of condensed notations concerns functional derivatives. 
We shall always reserve functional variational notations 5/5g a = 5/5g a (t) for varia- 
tional derivatives with respect to functions of time, while the variational derivatives 
with respect to functions of spatial coordinates will be denoted by partial derivatives 
d/dq l = 5/5g ab (^). 

By the same method we shall calculate g + -derivatives of the gravitational prefactor 
P and the higher-order derivatives of all the quantities in question. An important 
ingredient of all these calculations is the quantity 

dg a (t\q + ,q. 



It will be obtained from the following linear boundary value problem. 



(4.21) 



5. Linear boundary value problem 

To find ( |4.21| ) one has to consider the boundary value problem fl4.18| ) - ( (4.19| ) and 
make an arbitrary infinitesimal variation of the final point of the extremal 5q l + . This 
variation will induce the variation of the extremal 5q + ^ — ► 8g a = (5q l (t),SN li (t)) 
satisfying the linear boundary value problem 

8q(t+) = 5q + , 6q(t-) = . (5.2) 
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This follows from (4.18) after taking account that the varied trajectory must also 
be an extremum. The solution of this problem obviously gives the needed quantity 
(|4.21|) . However, the solution is not unique because of the local gauge invariance of the 
gravitational action ( 4.14j) . 

This gauge invariance with respect to spacetime diffeomorphisms, discussed above 
in the canonical formalism, in the Lagrangian formalism takes the form of the invari- 
ance with respect to infinitesimal transformations with local parameters = / M (£) 



and generators R a 



In view of the relation (|4.15|) between the canoni- 



cal and Lagrangian formalisms these gauge transformations can be obtained from the 
invariance transformations ( |3.10| ) - ( |3.11| ) [p4 ] 

(5.3) 
(5.4) 



g 
R) 



g a + r;f, 



S(U - t M ) 



dH„ 



dpi 



R„ 



R c , 



S 



d_ 



p=p°{q,q,N) 
X 



K s(u - 1 



p=dS/dq 



s(t a - tn). 



(5.5) 



Note that in the definition of gauge generators we used condensed notations in which 
they have the above form of delta-function type kernels with subscripts of time coordi- 
nates indicating to which spacetime condensed index they belong. Note also that gauge 
transformations of superspace coordinates are ultralocal in time in contrast with the 
Lagrange multipliers whose transformations include the time derivative of the gauge 
parameter, see (3.11). The invariance of action in terms of these generators takes the 
form of the following identities 



R" 



0. 



(5.6) 



M 5g a 

Their functional differentiation shows that on-shell, that is on the solution of classical 
equations of motion S a = 0, the Hessian of the action S ab is not an invertible operator 
because it has zero-eigenvalue eigenvectors - the gauge generators: 

6R a . 



R?,S, 



ab 



-s a - 



0. 



(5.7) 



This is of course a generic feature of constrained systems. The degeneracy implies that 
there is no unique solution to the linear problem ( |5.1| ) for Sg, which is always defined 
up to gauge transformations of the above type. The standard method of singling out a 
particular solution from their gauge equivalence class uses the gauge-fixing procedure 
[|T], This procedure in application to the full nonlinear equations of motion Ql.18] ) 
- ( [4.19|) consists in imposing gauge conditions 



x fl (g,g) = o, 



(5.8) 
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which should have a nondegenerate Faddeev-Popov operator (or operator of ghost 
fields) - the two-point kernel of their gauge transformation 

= J^aK, (5-9) 

det<X^0. (5.10) 



For the linearized problem ( |5.1| ) the gauge fixing can obviously be reduced to im- 
posing the linearized gauge condition on field disturbances Sg, 

Xt5g a = 0, X^jf- (5.H) 

In view of this condition, the linear equation ( |5.1| ) on 5g can be replaced by another 
equation 

F ab 5g b = 0, (5.12) 

F ab = S ab - X ^uXl (5-13) 

where the operator F ab is supplied by a gauge-breaking term with some invertible 
gauge-fixing matrix c^ u , cf. (3.21). In contrast with S ab , this operator is nondegenerate 
because now 

R»F ab = -Qfapxt (5-14) 

and the right-hand side represents here in view of (|5.10| ) a functional matrix with the 
rank coinciding with the dimensionality of the gauge group. 

Now we have to consider briefly the properties of the obtained operator F ab . Note 
that the original operator S ab was of second order in time derivatives only in the sector 
of g*-variables, Sik = 5 2 S/5q z (t)5q k (t'), because only these variables enter the gravi- 
tational Lagrangian with their velocities q. Correspondingly, the mixed sector = 
5 2 S/5q i 5N' J ' is a first-order differential operator in time, while S^ u = 5 2 S/5N fl 5N' / is 
ultralocal in time, that is proportional to an undifferentiated delta-function, S^ v ~ 
8(t[i — t v ). The addition of a gauge-breaking term in F ab can drastically change this 
structure, provided the chosen gauge conditions ( |5.8|) are so-called relativistic ones. 
Relativistic gauges contain time derivatives of lapse and shift functions so that, 

det^C ^ , (5.15) 

and the corresponding Faddeev-Popov operator is of second order in time derivatives. 
(This is a simple consequence of the fact that the gauge generator ( |5.5|) in the sector 
of Lagrange multipliers is itself a first-order differential operator.) It reads 

fi-yl 1 (ft 
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This choice of gauge conditions guarantees that the functional matrix Xa °f linearized 
gauge conditions is also a first-order differential operator with the delta-function type 
kernel 

xZ^id/dtMtp-ta). (5.17) 

In what follows we shall often denote the differential operators and the direction (right 
or left) in which they act by arrows. 

Thus, with the choice of relativistic gauge conditions the operator F a b is of second 
order in time derivatives acting in all sectors of test fields p> a = (<£> l , <£> M ). For any two 
such test fields (pi and (p 2 one can write down the Wronskian relations for the original 
Hessian operator S a b —Sab (d/dt)6(t — t') 



rf{t) Sab &{t) ~<f a i(t) Sab f\{t) 

d 



<Pl{t) Wab ~ V>l(t) Wab <f2(t) 



(5.18) 



and for the gauge-fixed operator F ab =F a b (d/dt)S(t — t') 

<pl{t)Fab fl(t) ~ <pt(t) Fab <f\(t) = 



d 

dt 



(5.19) 



In the Wronskian operator Wab (d/dt) only the components Wib are nonvanishing. 
They can be obtained from the gravitational Lagrangian by varying the canonical 
momentum conjugate to q l 

W ib (d/dt)5g\t) = -5 9L{q d ^ N) . (5.20) 

For the operator F^ all the components of its Wronskian operator are nonvanishing 
and related to Wab by a simple equation, 

Wab(d/dt) =Wab (d/dt) + -^c> X u b (d/dt). (5.21) 

Now we can construct the solution of the boundary value problem (|5.2| ), (|5.11|) and 
(|5.12 ) for a small disturbance Sg a of the classical extremal induced by a variation of its 
end point 5q\. First, let us introduce a special Green function D ab (t, t') of the operator 

F a b 

F ab D bc = 5 c a , (5.22) 

D ic (t±,t) = 0, (5.23) 

Xj£(d/dt±)D bc (t±,t) = 0. (5.24) 
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This Green function is uniquely defined by the Dirichlet boundary conditions on its 
g*-components ( |5.23 ) and by special Robin- type boundary conditions (|5.24 ). It plays 
the role of a graviton propagator analogous to the causal Feynman Green function, but 
with special boundary conditions appropriate to the definition of the two-point solution 
of the Wheeler-DeWitt equation ( 3.16|) . It was first introduced in [10] where it was 
shown that this is a propagator of the semiclassical loop expansion of this solution. 
Its boundary conditions have a number of remarkable properties including the BRST- 
invariance and selfadjointness of the graviton operator ||28|| . Here we add one more 



property which is proved in Appendix A: this Green function gives a solution to our 
linear problem for Sg a induced by 5q l + , 



6g a (t) = 8q\ W ih {d/dt + )D ba {t + ,t). 

From this follows immediately the desired answer for the derivative (|4.21 ) 
dg a (t) 



W ib D ba (t + ,t). 



(5.25) 



(5.26) 



The perturbation scheme can now be employed. 



6. Heisenberg operators of matter fields 



Now we can continue the calculation of the effective Hamiltonian ( 4.11|) in the 



modified Schrodinger equation (|4.6|) . Finding the derivative Q4.20Q requires the know- 
ledge of the functional derivatives of the unitary evolution operator Uq which is rather 
straightforward. In view of the chronologically ordered nature of the latter one has the 
expression (time arguments of evolution operators denote the time intervals in which 
they are acting) 

XT J numat 

" -iU (t + ,t)—-U (t,t_). (6.1) 



5g*(t) uv "' J dg«(t) 

This can be rewritten in terms of the Heisenberg operator of the matter Hamiltonian, 
HF*(t), as 

5U 



5g a (t) dg a (t) 

Here the Heisenberg operator H™ at (t) differs from the Hamiltonian H mat in the Schro- 
dinger picture 

Hr\t) = U Q {t+,t)H^U»\t+,t) = tf mat (y3(t),^(t),s(t)) (6.3) 
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by replacing the Schrodinger matter- field operators, (0,Pip), with the Heisenberg ones, 
(< ) 3(t),p v ,(t)). The Heisenberg operators satisfy operator equations of motion. There- 
fore, the derivative of the Hamiltonian in ( |6.2| ) with respect to the gravitational variable 
can be replaced by the functional derivative of the matter field action S mat (also taken 
in the Heisenberg picture)^ 

But the right-hand side of this relation is just the matter stress tensor in the Heisenberg 
picture of matter fields decomposed in the normal basis of (3+l)-foliation. Going back 
from condensed notations to the usual ones, it can easily be checked that 

<5S»« dg all A 1 „ 1 



T - = ra# = y 2 T {xl ^ TUxh -^™ x) 1 ■ (6 - 5) 



where 

X cmat 

T*{x) = 2—— (6.6) 

is a conventional matter stress tensor (density)Q. 
Thus we finally have 

' )Ut ' if a (t)U . (6.7) 



6g a (t) 

The substitution of this result into the g + -derivative of the evolution operator ( |4.20| ) 
gives rise to the quantity 

t a (t) = ~ f + &t' D ab (tX)f b {t') = -\D ab f b , (6.8) 
nip Jt- mp 

where we have again introduced a condensed notation in the last step. This quantity 
is a solution of the linear inhomogeneous equation (see Appendix A) 

rrip S ab i b (t)+f a (t) = 0, (6.9) 



5 Indeed, this follows from differentiating the relation L mat = p v ip — ff mat which gives extra terms 
proportional to the equations of motion for (p. Up to some operator-ordering ambiguity (usually 
absorbed by renormalization in renormalizable matter field theories), these terms vanish for Heisenberg 
operators satisfying these equations. 



6 The condensed notation ( |6 . 5| ) should not be mixed up with the notation (3.15) for the total 
momentum constraint in the canonical formalism. Since the latter will not be used below we take the 
liberty to reserve this notation for a spacetime covariant quantity - the stress tensor of matter fields. 
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which are the linearized nonvacuum Einstein equations with operator matter source. 
Therefore, i a {t) can be regarded as a gravitational potential generated by the back 
reaction of quantum matter on the gravitational background. It satisfies the same 
boundary conditions at t± as the Green function D ab and the same linearized gauge 
conditions 

xT = o. (6.10) 

As shown in Appendix A, these properties follow from the covariant conservation law 
for the Heisenberg operator of matter stress tensor ( |6.5[ ). This conservation is in turn 
a consequence of the diffeomorphism invariance of the matter action and Heisenberg 
equations of motion for matter fields 

r Cmat 

R% = = o (6.ii) 

(here R£ are the generators of local gauge transformations of matter fields - their local 
diffeomorphisms) . 

In terms of the constructed gravitational potential, the g + -derivative of the evolution 
operator ( 4.20p takes the form 

dtj — * 

= V m U = -i m 2 P W ma i a (t + )U Q . (6.12) 

The second-order (covariant) g + -derivative of Uq which is contained in the effective 
Hamiltonian Q4.11J) - the first correction term proportional to trip 2 - follows from 
differentiating (4.20) and using (6.7). It equals 

V m V n U = 1 \T> n — j T a Uo + J^E- (6-13) 

The second-order variational derivative of the evolution operator here can be obtained 
by repeated functional differentiation of ( |6.7|) to give the expression 



m " &-T(fjJf/o, (6.14) 



5g a 5g b 5g a 5g 

T (f a f b ) = 6{t a - t b )f a f b + 6{t b - t a )f b f a , (6.15) 

where T denotes the operator chronological ordering. 

Similarly, the first term in (|6.13| ) follows from repeated differentiation of (|5.26|) and 
expressing the answer in terms of functional variations of the gravitational background. 
In view of the variational equation for the Green function D ab 

SD ab = -D ac 5F cd D db = -D ac S cde D db 5g e , (6.16) 
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whereQ 

^ = (6 - 17) 
is a three-point gravitational vertex, one has 

{v m ^\ f a = -m%(V m W na ) i a + m 2 P (W ma D ac )(W nb D bd )S cde i e . (6.18) 



Combining ( |6.13| ) with ( |6.14| ) and (|6.18 ) and contracting these expressions with 



gmn one can g e ^. ^ e contribution of matter fields to the effective Hamiltonian (|4.11| ). 
It contains, however, the object {Q mn T> m Wna) i a which requires further simplification. 
The Wronskian operator here is a local quantity taken at the moment of time t + , 
but its g + -derivative calculated as above in terms of the background field variation, 
d/dq + = (dg a /dq + )S/Sg a , again produces a nonlocality - the Green function D ab (t,f) 
and its derivatives with respect to both arguments t and t' taken at the coincident 
points t = t' = t + . This can be condensely written down in the form of a special local 
three-point vertex w abc {t + ) 

{Q mn V m Wna) i a = D ab w abc (t + ) i c 

= D ab {t,t')w abc (d/dt,d/dt',d/dt'')i c {t'')\ t _ ti _ tn _ t . (6.19) 

The form of w a bc{t+) is rather complicated and can be found in Appendix B. 

With these notations the final form of the matter field contribution to - the 
first correction term - reads 

1 Q mn {V m V n U^)lJ^ = \m 2 P Q mn T ( W m J a W nb i b ) + )-D ab w abc (t + ) t c 



2mp v " ' ' 2 ^ V / 2 

~G mn {W ma D ac ){WnbD bd ){S cde i e + \ST) ■ (6-20) 
2 mp 

The resulting three terms can be given a Feynman diagrammatic representation 
with different structure. Note that because of (6.8) all terms are of the same order 
mp 2 . The first term begins with the tree- level structure quadratic in gravitational 
potential operators t a . Note that despite the fact that these operators are taken at one 
moment of time t + , their chronological product is nontrivial because it should read as 

T (t a i b ) = —j- D ac D bd T (f c f d ) (6.21) 



7 We assume for simplicity that the matrix of gauge conditions \a can b e chosen background field 
independent, whence SF cc i = S c d e Sg e - Otherwise extra vertices involving the variational derivatives 
SXa/^9 b wn l appear and generate an extra set of diagrams with ghost propagators of the Faddeev- 



Popov operator (5.9) 
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and, thus, includes all higher order chronological couplings between composite opera- 
tors of matter stress tensors. The second and the third terms in ( p.20 ) are essentially 
quantum, because their semiclassical expansions start with the one-loop diagrams con- 
sisting of one and two "graviton propagators" D ab , respectively. The quasi- local vertices 
of these diagrams are built from the Wronskian operators, gravitational three- vertices 
w abc{t+) an d S c d e and second-order variation of matter action with respect to gravita- 
tional variables 5^ at . The corresponding diagrams are shown in Fig. 1. 



7. Graviton loop effects 

The quantum gravitational contribution to the effective Hamiltonian ( |4.11|) is gen- 
erated by the gravitational preexponential factor P(q + ,qJ). In the one-loop approx- 
imation it is known as the solution fl3.23| ) of the continuity equation ( |3.18| ) which 
itself follows from the Wheeler-DeWitt equations. On the other hand, within the path 
integral representation of this solution one can write down the spacetime covariant 
representation of this quantity as a one-loop (gaussian) approximation for the path 
integral. Then, in contrast to the spatial functional determinants of (|3.23| ), it will be 



given by spacetime functional determinants of differential wave operators of gravita- 
tional and ghost fields (in a certain gauge- fixing procedure). This problem was solved 
T0|| , where it was shown that the prefactor is given by the one-loop effective action 



m 



of the theory calculated on the classical extremal fl4.18|) - ( |4.19| ) joining the superspace 
points q± at some moments t±, 

ir[g(t\q + ,q_)} 

T= ^TrlnF afe -iTrlnQ^ . (7.2) 

A remarkable property of this effective action T[g] is that its gauge (Fab) and ghost 
(Q$) inverse propagators coincide with the operators ( |5.13|) and (|5.9| ) introduced above. 



Moreover, their functional determinants are calculated with the same boundary con- 
ditions as those of the Green function D ab fl5.23j) - ([5.24 ) and Dirichlet boundary 
conditions of the ghost Green function Q~ lfl [[1(J. Therefore, to calculate the graviton 



contribution to Hf^ one can repeat the steps of the previous section and arrive at the 
equations similar to ( |6.12| ), 

V m P = -i m 2 P W ma i a (t + )P , (7.3) 



and flPOD 



' g mn {V m V n P)P- 1 = \m 2 p g mn (Wmal a )(Wnbl b ) + 1 -D ab w abc (t + )i c 



2m 2 p" y ~ '"' "~ '~ 2""*"' y " ' /v " ' ' ' 2 
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T cd = — — . (7.6) 



- l -Q mn (W ma D ac ) (W nb D bd )(S cde 7 e + \v cd ). (7.4) 
2 nip 

Comparing this with (6.20), the main modification is the replacement of the gravita- 
tional potential t a and the matter vertex by the new gravitational potential 7 a 
and the one-loop two-point "vertex" T cd , respectively: 

l a = -\D ab T b , r b = ^, (7.5) 
mp dg b 

5g c 5g° 

The new gravitational potential is generated by the one-loop stress tensor of vacuum 
gravitons r& which, similarly to ( |6.9D , enters as a matter source in the linearized Ein- 
stein equations for •y a . As is known [HJ (see Appendix A), this vacuum stress is also 
covariantly conserved on shell, 

R^T a = 0(5S/5g) (7.7) 

(which implies on-shell gauge invariance of the effective action), whence it follows that 
the new gravitational potential also satisfies the linearized gauge condition 

Xtr = 0. (7.8) 
Now the calculation of the cross term in H™ at presents no difficulty and gives 



1 g mn (V m P)p- 1 (V n U Q )U, 1 = m P Q mn (Wmal a )(W nb i b ). (7.9) 



The graphical representation of (7.4) is similar to Fig. 1, with i a (S™ d at ) replaced by 7" 
(T cd ), and no time ordering in the first graph. The graphical representation of (7.9) is 
very similar to the first graph in Fig. 1, with one t a replaced by r y a . 

8. Effective Hamiltonian and back reaction 



Collecting the equations fl6.20| ), ( |7.4| ) and ( |7.9| ) we get the total effective Hamilto- 
nian in the first order of the 1 / m^-expansion 

Hf = # mat + ^m 2 p g mn T (w ma h a Wn b h h ^ + l -D ab w abc (t + ) h c 

-\Q mn (WmaD ac ) (W nb D bd ) [S cde k + \{S™t + T cd ) ]. (8.1) 
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Here the full gravitational potential h a = i a + 7" determines the total back reaction of 
quantum matter and graviton vacuum polarization on the gravitational background, 

h a = — ^D^ifb + Tb). (8.2) 
mp 

It solves the linearized Einstein equations with the full stress tensor source and satisfies 
the gauge and boundary conditions of the above type, 

m 2 P S ab h b + f a + T a = 0, (8.3) 
X fh a = 0, fc'(tt) = 0. (8.4) 

This result was obtained in the first subleading order of the 1/mp-expansion, but the 
different terms in (|8.1|) have a very different nature from the viewpoint of a physically 
reasonable approximation scheme. As it has already been mentioned above, only the 
first two terms on the right-hand side have tree-level contributions, while the rest are 
essentially loop contributions. Formally this is explicitly indicated by an extra inverse 
power of mp. Certainly this is an artifact of the definition of the gravitational potential 
( |8.2| ) involving the inverse of mp and, therefore, formally all the terms belong to the 
same order of asymptotic expansion in Planck mass. However, from a physical point 
of view, one might consider quantum states with a big mean value of matter energy 
density, much higher than the energy of graviton vacuum polarization, so that 



5^ = 0(1), ^L = 0(l/m P ). (8.5) 
m P mp 

In this situation, then, only the back reaction from the matter part is relevant. This 
is, in fact, the situation most frequently studied in the quantum field theory on curved 
backgrounds. In this case of large matter sources, t a = 0(1) ^> 7°, only the first two 
terms of Hf remain dominating, and (|8.1|) reduces to 

Hf = # mat + \m 2 p Q mn T (w ma i a W nb i b ) + 0(l/m%). (8.6) 

Here we mainly restrict ourselves to this physical situation when the matter field back 
reaction is dominating over other effects. 

What is the interpretation of the second term in (8.6)? According to the definition 
of the Wronskian operator (|5.20|) , W ma i a is the linearized momentum conjugate to the 



gravitational potential t a , so that this term is just the kinetic energy of the gravitational 
radiation produced due to the back reaction of quantum matter sources. 

There is another important representation of this term which allows one to un- 
derstand better its dynamical properties and establish its relation to an analogous 
term found in H. This representation is based on a special decomposition of the back 
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reaction momentum W ma i a in the basis reflecting the gauge properties of the gravita- 
tional field. We know that the gauge direction in configuration space is defined by the 
generator ( gl2|) . By using the contravariant metric Q mn (|4.12|) and its covariant 
inverse, 

Qmn = ( S 1 )mn , (8.7) 

one can define in configuration space the orthogonal decomposition of tangent and 
cotangent vector spaces into subspaces longitudinal and transverse to the gauge gen- 
erator. The projector onto the transverse subspace looks as follows. Define the matrix 

iV> = V^ ifc V* (8.8) 

and assume that on the background of the classical extremal this matrix is nondegen- 
erate, denoting the inverse by 

jV^ = (jV-i)^. (8.9) 
Then the transverse projector equals 

IT™ = 8™ - V™A^V W , (8.10) 

v m = v™g mn . (8.ii) 

The distinguished role of this decomposition (with respect to such a metric) follows 
from the fact that the matrix flS.8| ) defines the matrix of second functional derivatives 
of the action with respect to Lagrange multipliers - lapse and shift functions 

= N u 5(t - f) (8.12) 

and turns out to be the coefficient of the variation of these functions in the linearized 
gravitational constraints H u (q,p°(q,q, N)). This property is used in Appendix C to 
show that the momentum of the gravitational potential Wmai a {t+) taken at the final 
moment of time t + has the following orthogonal decomposition 



Wmai a (t + ) = -Ii mn i n {t + ) - \ \7 mu N^f u (t + ), (8.13) 

nip 



where Il mn = Q mi Il l n . In view of this decomposition the effective Hamiltonian acquires 
the form 

Hf = H^ + i TO |n mn T(>K) +^ T N^f u (t + )f u (t + ) + 0(l/m 2 P ) (8.14) 

(this should not look confusing because T a (t + )/m 2 P = 0(1) and the third term is 
therefore dominating over 0(l/m p )). 



25 



An important difference between the second (kinetic) term and the third (potential) 
term concerns their locality properties. The kinetic term is essentially nonlocal, because 
in integral form it involves the fields at all moments of time t + > t > t_. On the 
contrary, the third (potential) term is ultralocal in all field variables taken at t + . 
Moreover, /^-components of matter stress tensor coincide (up to sign) with matter 
parts of the Hamiltonian and momentum constraints H™ &t {t + ). In addition, at t + 



the Heisenberg operators coincide with the Schrodinger ones (see the definition (|6.3|)). 
Therefore, this part of the effective Hamiltonian is just a quadratic combination, local 
in time, of the operators if° iat taken in the Schrodinger picture, 



2m P F 

This is exactly the structure in the effective iterational Hamiltonian captured in H0. 
This term quadratic in the superhamiltonian of matter fields was obtained there by 
calculating the "longitudinal" part of g + -derivatives of the above type. As we see, ex- 
tension to the case of the full configuration space and taking into account the "transver- 
sal" part of these derivatives results in the local quadratic term ( |3.15| ) plus the non- 
local kinetic term involving only the transversal velocities of gravitational potentials, 

It should be emphasized that the representation ( |8.14| ) and the canonical decom- 
position (^.13|) exist only for problems with the invertible operator (RTBJ) . This is 



an important restriction characteristic of the problems related to the thin-sandwich 
problem in which the nondegeneracy of ( |S.8|) guarantees the possibility to solve the 
constraint equations for fixed superspace coordinates q = 5w,(x) and their velocities 
q = 5w,(x). This property is violated, for example, for the linearized theory on flat- 
space background, or more generally on backgrounds with Killing symmetries [Q. For 
such backgrounds a similar decomposition should be modified by special techniques in 
the kernel of the operator ( |3.8| ) spanned by the Killing vectors of the background . 
This goes beyond the present paper and will be considered elsewhere. In contrast to 
the representation ( |8.14p , the original form of the effective Hamiltonian QB.6| ) remains 



valid on all possible backgrounds. 



9. Summary and Outlook 

The main result of our paper is the calculation of all quantum gravitational correc- 
tion terms to a given matter Hamiltonian up to order rrip 2 , the final expression being 

8 Cf. the first correction term in Eq. (42) in that paper. This term does not contain h explicitly, 
reflecting the tree-level nature of this term. The remaining, imaginary, correction terms in Eq. (42) 
of [0 are proportional to Ti and are contained in the imaginary terms of (8.1) in the present paper. 
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given in (8.1). This result both generalizes previous results |§ and gives an interpre- 
tation in terms of Feynman diagrams. Our discussion thus builds a bridge between 
the canonical and the covariant frameworks in the semiclassical approximation. (Since 
there is no spacetime at the fundamental level, this is all one can do in this respect.) 
In future applications we intend to apply these correction terms to concrete physical 
situation such as the one in JIB] , where a quantum gravity- induced energy shift to the 
trace anomaly in De Sitter space was calculated. 

There are, of course, still open issues. One important open point is to find a consis- 
tent regularization that preserves the ordering of |18|, |19| without producing anomalies. 
This is a contentious issue, and it was argued, for example, in [0 that anomalies nec- 
essarily occur which may even spoil the standard semiclassical approximation. 

The perturbative non-renormalizability of quantum general relativity can of course 
not be circumvented by the present approximation scheme. Thus, either full canoni- 
cal quantum gravity is consistent and the approximation scheme if taken at all orders 
becomes useless, or the full theory is inconsistent and one has to consider an alterna- 
tive approach such as superstring theory. But even in the latter case it might be very 
likely that the first quantum gravitational correction terms have the form presented 
in (8.1). This hope is based on the correspondence principle between the fundamental 
superstring theory and the low-energy classical gravity theory and local theory of renor- 
malizable gauge fields. The progress of these theories teaches that there should exist a 
sub-Planckian energy domain in which the predictions of covariantly regularized local 
quantum gravity consistently describe the low-energy limit of this fundamental theory 
of extended objects. Among the examples of applications in such an energy domain 
one can mention the recently proposed mechanism of generating the energy scale of 
inflation by loop effects in quantum cosmology BTI] - the phenomenon which, being on 



one hand essentially quantum gravitational one, on the other hand provides an effective 
suppression of sub-Planckian scales and thus justifies the semiclassical expansion. 

In our discussion we have considered states of the form (2.13), i.e. states where the 
dominant part is a phase obeying the Einstein-Hamilton- Jacobi equations. However, 
generally arbitrary superpositions of such states are expected to occur. In such a case 
there is no longer a unique background spacetime available as the starting point for 
the approximation scheme. In many realistic situations one can, however, understand 
how the various semiclassical components in the superposition become dynamically 
independent. The key mechanism is the process of decoherence by irrelevant degrees 
of freedom |32j . The results presented in our paper yield the necessary technical pre- 
requisites to study decoherence processes in the early universe at higher orders of the 
inverse Planck mass. It must also be mentioned that there might exist situations where 
already the lowest order of the semiclassical approximation breaks down [33[] . This will 
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be considered elsewhere. 

Finally, there is the question whether some of the quantum gravitational correction 
terms lead to an effective violation of unitarity in the matter sector. This is supported 
by the fact that imaginary terms occur explicitly, see (8.1). It has previously been 
shown that the occurrence of such terms would have important effects for the process of 
black hole evaporation However, it was argued that one can perform the splitting 
of the full wave function satisfying the Wheeler-DeWitt equation into gravitational and 
matter parts in such a way that this effective non-unitarity disappears |^5, 3f|. We 
hope to clarify this issue in a future publication. 



A Ward identities and properties of gravitational 
potentials 



We begin this Appendix by deriving the solution ( |5.25|) of the boundary value 
problem (|5.12 ), (|5.2j ). For this purpose we write first down the Wronskian relation 
(|5.19| ) with Lfl = 5g a (t) and <p\ = D bc (t,t'), integrated over t from £_ to t + . On the 
left-hand side we have Sg c (t') in view of the equations for Sg a (t) and D bc (t, t'), while the 
right-hand side gives only the contributions of upper and lower limits of integration: 



5g%t>) = ( Sg a W" ab D bc (t, t') - Sg a w" ab D bc (t, t> 



be I 



t=u 



t=t- 



(A.l) 



In view of the boundary conditions on the Green function Q5.23D - ( 5.24 ) and the 
relation (|5.21|) for the Wronskian operators of F ab and S ab , only one term 5q\ W% b 
D bc (t,t f ) survives on the right-hand side, because W \i b = 0, and 5g a by construction 
satisfies the linearized gauge conditions ( |5.11| ). This proves ( |5.25| ) and ( |5.26|) . 



Now let us prove that in ( |5.25| ) 5g a really satisfies the linearized gauge conditions. 
For this purpose let us use the Ward identity relating the gauge and ghost Green 
functions. Eq. ( |5.14[ ) implies that 



c a p X%(d/dt)D bc (t,t') = -Q- lp (t,t') R^d/dt'), 



(A.2) 



where Q~ 1/3 (t,t') is a ghost Green function subject to Dirichlet boundary conditions 
at t±, 



Q»Q* 1P = ^ Q- Q ^(t±,t') = o. (A.3) 

Now acting on (|5.25|) by \a an d using (|A.2 ), one can see that in view of the differ- 



ential structure of the gauge generators ( |5.4j ) and ( |5.5| ) and the Wronskian operators, 
the time derivatives acting on the ghost Green function at t + cancel and the rest of 
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the terms vanish due to the Dirichlet boundary conditions for Q a lf3 , whence it follows 
that 

dq a 

X-| = 0. (A.4) 



The linearized gauge condition (|5.11 ) is thus fulfilled. In proving this we used the 



expressions for Wronskian operators 

Wik (d/dt) = Sik-^ + local terms, (A. 5) 



7 &L(q,q,N) dH^(q,p°(q, q, N)) , 

W ^ = ~ dq*dN» ~ ft? " GlkV » ■ (A ' 6) 



Let us now turn to the gauge property of the gravitational potentials t a and 7 a . 
Acting on (p.8|) by the gauge matrix and using the Ward identity of the above type, 



one can see that the result is proportional to R^Ta which is zero in view of ( 6.11 ). This 
proves (|6.10|). Similar proof holds for the gravitational potential r y a . 



B The vertex w a b c (t + ) 

The derivation of the vertex begins by noting that in view of (|A.5 ) and ( A.6|) 



— > dt n (t ) 

Wmat a (t + ) = -g mn ) +! + V m J»{t + ) . (B.l) 

dt± 



Therefore, the calculation of the derivative on the right-hand side of ( |6.19|) in terms of 



the background field variation, d/dq + = (dg a /dq + )S/Sg a , takes the form 

(G mn V m Wma) t a = - Wma D a± (t + , t+) GT Wnb t\t + ) 

+ (5I?f/5g b ) D bc (t + ,t + ) Wcm r(t + ), (B.2) 

where the variation of V mfl in ( |B.1| ) is replaced by the variation of the Lagrangian 
generator R™ = V™| p=p o, and it was taken into account that the covariant superspace 
derivative of the three-entry object G™ n is zero. The kernel of the first-order differential 

operator (5R™/5g b ) is understood here as acting on the first argument of the Green 
function. Thus this expression really reduces to the one-loop tadpole structure - the 
coincidence limit of the graviton Green function and its derivatives - and it can be 



represented by a quasi-local vertex w a b c (t + ) read off (|B.2|) . 
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C Gauge decomposition of gravitational potentials 



The proof of the gauge decomposition of the momentum conjugated to the gravita- 
tional potential ( |3.13| ) is based on (|B.!| ) above. The /i-component of the gravitational 
potential fl6.8| ) there involves the D fia (t + , t) -component of the Green function. It can 
be found from the /^-component of Q5.22| ) 



S»iD ia (t+,t)+ S,„D» a (t + ,t) =X«c aP X%D ba (t + ,t). (C.l) 

Let us use the Ward identity (|A.2| ) and take into account the Dirichlet boundary- 
condition on the ghost Green function. The latter implies that on the right-hand side 
here only the time derivative part of X ® = — (dx a /9N IJ ')d/dt + . . . will survive acting 
on the ghost propagator at t + . Similarly, on the left-hand side the operator coefficients 
Sfj,i and Snu are the first-order differential operator and ultralocal operator given by 

-* d 2 L(q,q,N) d d . . . 

Sfli= dtfdNi* dl + localterms = ~ V ^J t + localterms ( c -2) 

and (p.!2| ), respectively. According to our assumption, the functional matrix N^ u in 
( |8.12| ) is invertible, so that the equation resulting from ( |C.1| ) can be solved with respect 
to D va {t+,t), 



D-"(t + ,t) = AT«"V„ m — D""(t+,t) + N""^L.— Q- a v, (U,t) R %dldt). (C.3) 



Note that the a = a-component of the generator Rp (Rp —d/dt 5p + . . .) is acting 
to the left in this equation. When contracting it with T a to obtain one should 

integrate in the last term by parts in order to use the conservation of matter stress 
R pT a = which yields 



dt m 1 f)-v a d 

f(t + ) = N""V™ N i»™«Q-W {t t)T ( t) 

dt + m P dN u dt + 



(C.4) 

t=t + 

The contribution of the surface term at t + here is not zero even for the Dirichlet 
boundary conditions, because the coincidence limit at t + is understood in the sense 
when the second argument of (d/dt + )Q~ a 1/3 (t + ,t / ) tends to t + after the first one. It is 
easy to show that 
dx a d 

- (/ . i T + . I ] 

t=t+ 



dN» dt. 



Ql 10 (t + ,t) _ = W , a {d/dt + )Q- a ^{t + ,t) 



t=t 4 



(C.5) 



where W ^ a (d/dt + ) is a Wronskian operator corresponding to the Faddeev-Popov op- 
erator ( |5.9| ). The second equality here can be easily proved from the spectral decom- 
position of the ghost Green function subject to Dirichlet boundary conditionsP] p9|. 



9 Note that this property is a direct analogue of the following equality for the matrix J5.26| ): 
dq k (t+)/d<t + = 8*. 
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Thus 



Aim U \ i 

t»( t+ ) = N^V um —^- -—N^t(t + ). (C.6) 



Substituting this relation into (BT) we finally get the decomposition ( 8.13| ). It is obvi- 



ous that the same proof holds for the gauge decomposition of the one-loop gravitational 
potential r y a . 
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Figure Caption 

Fig 1. Contributions of quantum matter fields to Hf s , see (6.20). The time parameter 
t+ labels the vertices at the spacetime point at which Hf 1 is evaluated. Dashed lines 
labelled by T denote the chronological ordering of matter stress tensors in the bilinear 
combinations of gravitational potentials t a . 
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